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0^ ■ We refine Misner's analysis of the classical and quantum Mixmaster in the fully inhomo- 

■ geneous picture; we both connect the quantum behavior to the ensemble representation, 

both describe the precise effect of the boundary conditions on the structure of the quan- 
ta "3 ' turn states. 

Near the cosmological singularity, the dynamics of a generic inhomogeneous cos- 
mological model is reduced by an ADM procedure to the evolution of the physical 
degrees of freedom, i.e. the anisotropics of the Universe. In fact asymptotically to 
the Big-Bang the space points dynamically decouple 1 because the spatial gradients 
(3Jq| of the dynamical variables become of higher order 2 and we can model the inhomo- 

geneous Mixmaster via the reduced action: 



I = J d 3 xd,T (p u d T u + p v d T v - e) , t = v\/p 2 u+pl, (1) 

The dynamics of such a model is equivalent to (the one of) a billiard-ball on a 
Lobachevsky plane; this can be shown by the use of the Jacobi metric. 1 The man- 
ifold described turns out to have a constant negative curvature, where the Ricci 
scalar is given by R = —2/E 2 : the complex dynamics of the generic inhomogeneous 
model results in a collection of decoupled dynamical systems, one for each point of 
the space, and all of them equivalent to a billiard problem on a Lobachevsky plane. 
We want to investigate the relation existing between the classical and the semiclas- 
sical dynamics, and from this analysis we will derive the correct operator ordering 
to be used when quantizing the system. 

Let's write down the Hamilton- Jacobi equation for the system 



(2) 
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This can be explicitly solved by separation of constants 



Sq(u, v) = k{y a )u + ^e 2 - k 2 (y a )v 2 - eln 2 * = + c(y a ) . (3) 




The semiclassical analysis can be developed furthermore, and the stationary con- 
tinuity equation for the distribution function can be worked in order to obtain 
informations about the statistical properties of the model; as soon as we restrict 
the dynamics to the configuration space, we get the following equation for the dis- 
tribution function w 



dw(u,v;k) U E\ dw(u,v;k) E 2 — 2k 2 v 2 w(u,v;k) , 
+ VVW dv + fcw 5 y/E* - (kv) 2 = U 

This can be solved, and the exact distribution function can be obtained as soon as 
we eliminate by integration the constant k 



1 



w(u, v) = / — — — —dk (5) 

J -EL vvE 2 — k 2 v 2 

It is worth nothing how in the case g — const, the microcanonical Liouville measure 
w mc (u,v) = is recovered. 

We expect that the distribution function w(u, v) is re-obtained as soon as the quan- 
tum dynamics is investigated to the first order in h. This can be easily done as a 
WKB approximation to the quantum dynamics is constructed; as soon as we retain 
only the lower order in h, we obtain that: i) the phase S(u,v) coincides with the 
Hamilton-Jacobi function, and ii) the probability density function r(u, v) obeys the 
following equation: 

k 9 -^ U l(*) 2 -k 2d -^l + ^ - * V) " ^ r(u, „) = . (6) 
du V \ v J dv v 2 VE 2 -k 2 v 2 

that coincides with ([¥]) for a particular choice for the operator-ordering only, i.e. 

With this result, the problem of the full quantization of the system can be taken 
in consideration. The main problem is the presence of the root square function 
in the definition of the Hamiltonian |T]), but well grounded motivations exist 3 to 
assume that the real Hamiltonian and the squared one have same eigenf unctions 
and squared eigenvalues. This way the solution of the eigenvalue equation can be 
obtained 

^(u, v) = a n K s _i/2{^nirv) sm(2nitu) (8) 

n>0 

The spectrum is obtained as Dirichelet boundary conditions on the domain Tq (see 
FigHJ) are taken into account. The condition on the vertical lines can be imposed 
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exactly, but the one on the semicircle cannot be solved exactly; so we approximate 
it as in FigOwith a straight line v = 1/ir (it preserves the domain area /i = tt). 
All these imply that s = 1/2 + i, with t £ 5i, and that the spectrum assumes the 



7 . I 

Fig. 1. The real domain Fq Fig. 2. The approximate domain 

following form 

Et = ( 1 -+t 2 )H 2 (9) 

The values of the real parameter t have to be numerically evaluated by solving the 
equation Ki t (2n) = for every natural n. This condition implies a discrete but 
quite complicated shape for the spectrum. 

Asymptotic expansions for high occupation numbers can be derived for different 
regions of the parameters (t, n) A 
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